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We study the general properties of the CMB temperature four-point function, specifically its 
harmonic analogue the angular trispectrum, and illustrate its utility in finding optimal quadratic 
statistics through the weak gravitational lensing effect. We determine the general form of the trispec- 
trum, under the assumptions of rotational, permutation, and parity invariance, its estimators on the 
sky, and their Gaussian noise properties. The signal-to-noise in the trispectrum can be highly con- 
figuration dependent and any quadratic statistic used to compress the information to a manageable 
two-point level must be carefully chosen. Through a systematic study, we determine that for the case 
of lensing a specific statistic, the divergence of a filtered temperature- weighted temperature-gradient 
map, contains the maximal signal-to-noise and reduces the variance of estimates of the large-scale 
convergence power spectrum by over an order of magnitude over previous gradient-gradient tech- 
niques. The total signal-to-noise for lensing with the Planck satellite is of order 60 for a ACDM 
cosmology. 
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I. INTRODUCTION 



The power spectra or two-point correlations of Cosmic Microwave Background (CMB) temperature and polarization 
\ anisotropies are widely recognized as a gold mine of information on cosmology. These spectra in fact contain all of 
J> ■ the information embedded in the CMB if the underlying fluctuations are Gaussian distributed. Nonetheless even if 
C \ the initial density fluctuations are Gaussian, non-Gaussianity in the CMB temperature fluctuations will be generated 
by non-linear processes. These generally are associated with the secondary anisotropies that are imprinted as the 
\ photons propagate through the large-scale structure of the Universe from the epoch of recombination. 

Secondary signatures in the three-point correlation of temperature anisotropies have recently received much at- 
y—i \ tention [1-3] following early pioneering work on intrinsic correlations in the initial conditions [4,5]. The four-point 
• correlation and its harmonic analogue the trispectrum has received considerably less attention despite the fact that 
it directly controls the noise properties of estimators of the power spectrum. In particular, an all-sky treatment of 
O ,■ the trispectrum that incorporates the full rotational symmetry properties of the trispectrum has been lacking in the 
literature (c.f. [6]). Exploitation of the symmetry properties can assist in the isolation of the physical mechanisms 
; I ■ underlying the generation of the trispectrum as we shall see. 

In this paper, we establish the framework needed to study the trispectrum on the full sky. We begin in §11 with a 
discussion of the symmetry properties of the n-point function on the sky, with an emphasis on the 4-point function, 
| and their implications for the general form of the harmonic spectra. We consider estimators of the trispectrum and 
their noise properties in §111, and the trispectrum-based power spectra of quadratic statistics in §IV. Calculational 
techniques and relationships to the flat-sky formalism are given in two Appendices. In §V, we consider the specific 
case of the trispectrum generated by weak gravitational lensing of CMB photons by the large-scale structure of the 
Universe and show that there exists a quadratic statistic that optimally recovers the projected gravitational potential 
power spectrum (or convergence) on large scales. We conclude §VI. 



II. SYMMETRIES 



In this section, we derive the requirements that rotational, permutation, and parity symmetry impose on the 
n = (2, 3, 4)-point correlation functions on the sphere and their spherical harmonic analogues: the power spectrum, 
bispectrum and trispectrum. We begin with general considerations for the n-point function in §11 A, review the 
implications for the power spectrum and bispectrum in §11 B, and derive the consequences for the trispectrum in 
§11 C. In §11 D, we show how to construct trispectra with the required symmetry properties. 



A. General Considerations 



We begin by requiring statistical isotropy of the n-point correlation function on the sphere and its harmonic analogue 
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(e(n 1 )...e(n„)}= E E (e Iimi ...e lBTOB )y I - i (n 1 )...i-'-(n n ). (1) 

li...l n mi...m n 

Statistical isotropy demands that the n-point function is invariant under an arbitrary rotation R whose action on a 
spherical harmonic is expressed in terms of the Wigner-D function 

R[Yr (n)] = J2 Dl ™>m («, 0, (n) , (2) 

m' 

where a, (3 and 7 are the Euler angles of the rotation. To obey rotational invariance the harmonics must obey the 
relation 

(e, iroi ...e, Bron )= E (e lim[ ...e lnm/ jD^ im ,...D l ^ m , n , (3) 

for all a, /3, and 7. The reduction of this relation proceeds as follows. Each pair of rotation matrices may be coupled 
into a single rotation via the group multiplication property (or equivalently the addition of angular momentum) 

/ l\ I2 L \ f l\ I2 L 



D t im [ D t 2 m' 2 ~ E ( m \ m 2 2 - M ) ( m\ m' 2 —M' ) ( 2 ^ + 1)(-1) M+M D L MM , . (4) 

LMM' V / \ 1 -i / 



When the product is reduced a pair of D matrices, one seeks the form of the harmonic n-point function that reduces 
the pair to the orthogonality condition for rotations 

E(_i \m 2 -m pfh p)h _r (r\ 

m 

which is valid for an arbitrary rotation. The indices can then be permuted to find alternate orderings of the pairings. 
Invariance under a parity transformation which takes n — ► n 

vr - {-i)% m , (6) 

would require that 

n 

E h — even. (7) 

Reality of the underlying 6 field and the fact that 

Y ™* = (-l) m Y- m , (8) 

requires that 

©r* = (-i) m er m . (9) 

B. Power Spectrum and Bispectrum 

For the 2-point function there is only one step. The reduction of equation (5) requires the form 

{ < dl 1 m 1 Ql 2 m 2 ) — ^i/2^mi-m 2 ( — l) mi C/i ■ (10) 

For the 3-point function one first collapses one product of rotation matrices leaving 

(e Iiroi ...e l8m8 )= E ( e 'i-i ••• e W E ( 2i + i)(-i) M+M ' (ii) 

m' 1 ...m' 3 LMM' 

( l\ I2 L \ I li I2 L 



V mi m 2 —M J V m[ m' 2 —M' ' D MM'^t. m >. ■ 
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In order to reduce this relation to the orthogonality condition Eqn. (5), the sum over m! of the three-point function 
must be proportional to Sli 3 5m' -m' ■ Recalling the identity 

E(h h L \ / l\ l 2 L' \ _ 8ll>5mm' 
1 mi m 2 M / I mi m 2 M' ) ~ 2L + 1 ' ^ ' 

we can obtain the desired relation if the m dependence of the 3-point function obeys 

<e Iiroi . . . e hm3 ) = ( ^ ^ ^ J B hl2h . (13) 



C. Trispectrum 

The form of the 4-point function follows the same steps except that we use the group multiplication properties to 
pair say (li,h) and (l3,h) leading to the condition 

L 3i M 3i M' 3i V / \ 4 J4/ 

The same reasoning that led to the choice of the form of the three-point function implies that the following form is a 
solution 

< e - e <-> - E (A i -m) (i i m) <-i>"«B:m- («) 

Geometrically, Qi fJL) represents a quadrilateral composed of sides with length l\ . . . I4. The index L represents one 
of the diagonals of the quadrilateral and is also the shared third side of the two triangles formed by the corresponding 
pairs of sides. The Wigner-3j symbols in Eqn. (15) ensure that the triangle inequalities are satisfied. For this reason 
we will often refer to a set h, h, h, L) as a given "configuration" of the quadrilateral. 

The two other unique pairings of the indices, {h,h) and (h,h), yield alternate representations of the 4-point 
function. These are not independent since all three couplings yield complete sets according to the theory of the 
addition of angular momenta. The alternate representations are constructed as linear combinations of the (h,fo) 
representation with weights given by the Wigner-6j recoupling coefficients (see Appendix) 

Q^W = E(-i)^(2L + i){^ I L l}q\£(l>), 

Q^(£)=£(-l) i+i '(2£+l){ h ^ } Ql ^ L '^ ^ 

where we have used Eqn. (A4) to project one coupling scheme onto another. 

Symmetry with respect to the 4!/3 = 8 remaining permutations (2 orderings of the pairs, 4 ordcrings within the 
pairs) requires that 

QlllKL) = (-lf"Q\l\l(L) = (-l)^Q\±{L) = Q\f 2 {L) , (17) 

where Ej/ = l\ + I2 + L and = Z3 + I4 + L. If the four-point function is parity invariant then 

Q l i±(L)=Q l l ±(L). (18) 

We shall show how to construct trispectra that obey these properties in the next section. 

Finally it is useful to separate the contributions from the unconnected or Gaussian piece and the connected or 
trispectrum piece 

Q\i 2 SL) = G\±{L)+T%l{L), (19) 

where 

G\fc(L) = (-l) (l +V(2*i + 1)(2Z 3 + l)C h C h S hl2 S hU S L0 + (2L+l)C h Q 2 [(-l) l ' +, ' +L 5 h i a 5 hU + 8 llU 8 hh ] . (20) 
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D. Enforcing Symmetries 



The symmetries of the trispectrum described above may be enforced by the following construction. First describe 
the four-point function by a form that is explicitly symmetric in the three unique pairings 

(e iimi ...e, ro4 ) c = ^^(L)(^ £ _ L M ) £ ^)(-D- + (^/s) + (^/ 4 ), (2i) 

where c denotes the fact that we have removed the Gaussian piece of Eqn. (20). The two latter pairings can be 
projected onto the (I1J2) basis with the help of the Wigner-6j symbol to give 

Tl£(L) = P}£(L) + (2L + l)Y^ 
Within the three unique pairings, there are 4 permutations of the ordering implying that P is constructed as 

= Kt + (-v^t + (-v^Kt + (- i ) Sc,+Sl ^1 i • ( 23 ) 

The reduced function T underlying the trispectrum is an arbitrary function of its arguments except that it must be 
symmetric against exchange of its upper and lower indices 

Kt( L )=Kt( L )> ( 24 ) 

and if parity invariant obeys 

T^{L)=T^{L). (25) 
This then completes the enforcing of the rotation, permutation and parity symmetries of the trispectrum. 

III. ESTIMATORS AND SIGNAL-TO-NOISE 

We show in §111 A that the fundamental estimator of the trispectrum involves a weighted sum over the multipole 
moments in a given quadruplet of harmonics. These estimators have well defined noise properties as derived in §111 B 
which can be used to calculate the theoretical signal-to-noise in the trispectrum. A non-vanishing trispectrum can on 
the other hand decrease the signal-to-noise in the power spectrum by introducing a covariance between its estimators 
as shown in §111 C. 

A. Estimators 

From the orthogonality properties of the Wigner-3j symbol, one can invert the relationship for the four-point 
spectrum to form the estimator 

iffiW-pi + i) Y. H)"(i^)(: , J i_ I M )(e„.,..e ta )-fii;!:, (26) 

where the estimator for the Gaussian piece is constructed out of those for the power spectrum. Note that for 
configurations whose I sides are not equal in pairs, the Gaussian piece vanishes and the sum over m's of the spherical 
harmonic coefficients is an unbiased estimator of the trispectrum. 

We can alternately form an estimator of particular configurations of the trispectrum directly from the sky map 
itself without an explicit expansion in spherical harmonics. Following Spergel & Goldberg [8], let us define a new set 
of sky maps weighted in rings centered around a point q: 

e«(q) - / rfA0 (") P <(" • • (27) 

Expanding the Wigner-3j symbols in terms of spherical harmonics and using the addition theorem, we obtain 
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(22) 
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f%l{L)+G\±{L) 



(2L + 1) 



f dCjg r 

J 4tt J 



dqb 

47T 



eii(qa)e; 2 (q a )ei 3 (q b )e (4 (q b )P L (q a • q b ) . 



(28) 



Since the Wigner-3j symbol vanishes if l\ + I2 + L = odd, this expression can only be used to estimate even terms. 

To measure all configurations of the trispectrum is, needless to say, a daunting task. Aside from the computational 
expense, one must also treat complications associated with estimators of harmonics on a fraction of the sky. Even for 
an all-sky CMB experiment, the removal of galactic foregrounds will limit the data to a smaller fraction of the sky 

/sky ■ 



B. Signal-to-Noise 

Returning to the estimator of Eqn. (26), one can calculate the Gaussian noise variance of the estimator, 

(fl£*(L)fl£(L')) = (2L + l^Cjf CjfCjf Cjf , (29) 

if no two Vs are equal. Here Cf ot is the sum of all contributions to the power spectrum including the intrinsic CMB 
fluctuations, instrumental noise and residual foreground contamination. 

From the permutation properties of Q (or T) in Eqn. (16), the full covariance of the estimators then becomes 



(f}£*{L)0l?(L')) 



= 5 LL ,5\l + (2L' + 1) 



1 1 \h+i 3 \ n h L 1 13 



f h 


ro 


L 1 r!4 




U 


V ) 632 



(30) 



(2L + ljCf^Cf^Cf^Cff 

if no two I's in the primed and unprimed sets are equal. Here 

C = [Shiji 2h + (-lf u 5i lh Si 2 i a ][6 h iJ l4 i d + (-lf L S hld S ulc ] + [a <-> c][b <-> d] , (31) 

accounts for the permutations within the three fundamental pairings. Recall that = I1 + I2 + L and = + + L 
The two terms involving the Wigner-6j symbol reflect the fact that alternate pairings of the indices supply redundant 
information in both the signal and the noise. 

If any two Vs are equal, then the covariance has extra terms associated with the internal pairings in the primed 
and unprimed sets. Based on the fundamental relation 

(^(Ljf «(L')> = (-l)' 1+r ^Lo^o yj(.2k + l)(2/i + l)[5 l3 v 3 5 liVi + (-l)^S l3l ,5 hll JC h C,C l3 C l4 , (32) 

other pairings can be found through the permutation properties of Q (or T). No fundamentally new terms are 
introduced if three or four Z's are equal but each set of possible internal pairings in the primed and unprimed sets 
must be separately accounted for. 

The total signal-to-noise for each L in the four-point spectrum is 



(|) 2 ^ E E (Tti:%L))\Cov-'}(T$(L<)) 

^ ' hhhULl'l'l'l'L' 



2L + 1 C, tot C, tot C, tot C, tot ' 1 ' 

l h>i 2 >i 3 >U h h h u 

where "Cov -1 " indicates the matrix inverse, with elements labelled by their configuration (h, I2, h, U, L), of the 
covariance in Eqn. (30). In the second line, the restricted sum eliminates the 4!=24 redundant permutations above 
and neglects the signal-to-noise contributed when the the Vs arc equal. In the high signal-to-noise regime, one must 
also include the sample variance of the signal. On a cut-sky, the considerations of Appendix B imply that the overall 
(S/N) 2 is reduced by a factor of /sky- 
Note that if the sum in Eqn. (33) is not restricted the covariance supplied by the alternate pair orderings in Eqn. (16) 
necessarily contains off diagonal terms that mix L and L' . The covariance is distributed across many L's and can 
lead to overestimates of the signal-to-noise in 4-point related statistics by a factor of y/3. 
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C. Power Spectrum Covariance 



The trispectrum can affect two-point or power spectrum statistics by introducing a covariance between the estima- 
tors. The covariance of power spectrum estimators averaged over m is given by 

[Cov] ili2 = — - - — - ^ E ( iim 1 6r imi e '2m 2 6r 2 m 2 > ~ (^1X^2) 



1 



^(-l) h+h Q\±(0)-(C h )(C h ). 



^2h + 1 V2?2 + 1 

The expression for the covariance can be further broken into its Gaussian and non-Gaussian pieces 
2 



(34) 



[Cov] W2 = — — C h Si lh 



V(2Zi + 1)(2Z 2 + 1) 



h+h 



V(2/i + l)(2/ 2 + l) 



E {-l) L Pl±{L) 



L = \h-h\ 



(35) 



The effect of covariance for the signal-to-noise for the estimation of a set of underlying cosmological parameters pi 
can be calculated through the Fisher matrix 



\^ dC h ! dQ 2 



(36) 



where "Cov x " indicates the matrix inverse of the covariance in Eqn. (35). In particular, if the only parameter of 
interest is the overall amplitude A of a known template shape, then Faa — (S/N) 2 (see [3]). 



IV. POWER SPECTRA OF QUADRATIC STATISTICS 



Measuring all of the configurations of the trispectrum or four-point function is a daunting challenge. In this section 
we consider statistics based on the identification of points in pairs in the four-point function. These quadratic statistics 
may be optimized in signal-to-noise for their power spectra by filtering the original temperature field. We begin with 
general definitions for the quadratic fields in harmonic space (§11 A) and continue through a discussion of filters (§IVB) 
to a consideration of specific quadratic statistics (§IVC-IVH) and related cubic statistics (§IVI). The specific statistic 
and filter set that optimizes the signal-to-noise will depend on the configuration dependence of the trispectrum signal 
that is to be extracted. 



A. General Definitions 



To probe various aspects of the trispectrum, we can form the two-point or power spectrum statistics of a quadratic 
combination of the underlying field. To enhance the signal-to-noise, we begin by filtering the fields before collapsing 
the configuration, 

Q a (n) = J2e lm ftYr(n), (37) 

Ira 

where the index a = 1,4 to allow for 4 independent filters on the fields. In general, the identification of points in 
pairs implies that each pair (ab) involves a quadratic combination of the filtered field which in turn involves a mode 
coupling sum of the harmonic coefficients 

< M - (-D M E E <AL)e limi e l2m2] j^ £ L M ) , (38) 

hm.2 

where 



G 



(L) = /£ fl V(2*i + l)(2Z 2 + lK« 2 (i) , (39) 

and xi ± i 2 (L) represents different weights for different statistics x as specified below. The power spectra statistics 
relating two general quadratic statistics x and x may be separated into the non-Gaussian signal and Gaussian noise 
as 

<4m^V> - 5ll>5mm> (Cf + Nf) , (40) 

where 

Cf ^^2rri £ 4 2 r 2 (i)^3l( i )(- l ) ii+ ' 2+L ^t( i )' ( 4l ) 

and the Gaussian noise is 

N L* = ^H x hh{L) [xfX(L) + {-l) l ^ +L x^ h {L)] CjfCjf . (42) 

hl-2 

It will be useful in the following discussion of noise variance to also define the following two auxiliary power spectra, 

v l h12) = ^ E W W + (-i)' 1+i2+i 4l W] cr eg* . 

il(2 

^(34) = ^_ £ (L) ^34^ (£) + { _ l)h+l2+L ^ h {L)] Q to t Q tot (43) 

The signal-to-noise ratio in this power spectrum statistic can be calculated from equation (30) for the covariance 
of the trispectrum, 

(Cf Ct*') « ^— [ { y£*'< 12 >><vf' (34) ) + (Nf){Nf)] , (44) 



such that 



S\\ y. (CfYfiTX_ (45) 

Z—/ iQXXQx'x'\ ' ^ ' 

Lxxx'x' \ L LI 



N 



Strictly speaking, this is an inequality since we have neglected the covariance between Us dictated by the trispectrum 
covariance Eqn. (30). Since the trispectrum covariance is distributed broadly in the allowed L's, this signal-to-noise 
estimate is reasonable if we restrict the range of interest in L to a small fraction of the allowed range. The covariance 
can at most reduce the total signal-to-noise by a factor of \fi for the 3 unique pairings in the trispectrum. 

B. Filters 

The filters and specific form of the statistic x can be chosen to eliminate Gaussian noise bias and/or maximize 
the signal-to-noise. If (fl, ff) and (/f, ff) do not overlap in I, then the Gaussian noise bias of Eqn. (42) vanishes. 
Furthermore, trispectrum covariance between differing L's is identically zero and Eqn. (45) becomes a strict equality. 

For example the two filters may be band limited in mutually exclusive bands or parity limited 

e e (n) = i[0(n) + e(-n)] , 

f i ,2 _ f 1, I =even, , . 

/ ' =/ '=\0, Z=odd, (46) 

9 (n) = \ [6(n) - 0(-n)] , 

,3 _ fi _ J 0, I =even, , . 

h -Ji=\ 1 i=odd. (47) 
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This choice does not eliminate the auxiliary variance power spectra in Eqn. (43) and more generally does not maximize 
the total signal-to-noisc. Only if the filters are equal in pairs j\ = ff and ff = ff are the noise and auxiliary variance 
power spectra equal such that 

(CfCr') « {(Nf)(Nf) + (Nf)(Nt)) , (48) 

becomes the familiar form for the variance of the power spectra of a set of Gaussian random fields x. 

A comparison of the signal-to-noise in the full trispectrum Eqn. (33) and in a particular quadratic statistic x 
Eqn. (45) shows that the latter approaches the former if 



hh 

<iX(L)x% u (L)^w(L) gtotc tt cto ; c tot . (-W ) 
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where w(L) is an arbitrary function of L. To the extent that the right hand side is factorable in l a , a — 1, 4 the filter 
functions /" can by chosen construct this optimal statistic. Since the trispectrum is in general not factorable, we will 
next consider a wide range of choices for the quadratic ^-statistic which can be used to construct optimal statistics 
for various types of trispectrum signals. 



C. Temperature- Temperature 

The simplest quadratic statistic that we can form is the product of the filtered temperature field itself, 

e a (n)6 h (n) ee s ab (n) = £ sf M Y™ "(n) , (50) 

LM 

where s L b M is given by the general prescription of Eqn. (38) with x = s and the weighting 

s h i 2 (L) = (o o o) ' 6Ven ' ^ 

"Even" denotes the fact that s selects out li + I2 + L =even by virtue of the Wigner-3j symbol. The non-Gaussian 
power spectrum C£ s is then given by Eqn. (41) in terms of the trispectrum. The total signal-to-noise of this statistic 
can be estimated by retaining just the x = x' = x = x' = s terms in Eqn. (45). 



D. Temperature-Gradient 

The product of the filtered temperature field and the gradient of the filtered temperature field probes another aspect 
the trispectrum. This product is a vector field on the sky and may be broken up into components as 

e^V^n) ee ]T -j=[ ai ± *« 2 r b (n)i=(e T e e ) t . (52) 

The components a\ ± icti are spin-1 objects that can be decomposed in the spin-1 spherical harmonics [9], 

[a, ± za 2 ] ab (n) = ^(c ± ig)f M ±1 r L M (n) , (53) 

LM 

where c and g are the multipole analogues of the curl and gradient pieces. These quadratic statistics again follow the 
general form of Eqn. (38) with x = c, g and weightings 

c hl2 {L) = -^yJ~Tl)^ \ , odd, (54) 

9 h l 2 (L) ee iy/h(h + l)^ \ J ) , even. (55) 

Here and below "even" ("odd") denotes the fact that the expression holds for l\ + 1 2 + L =cven (odd) and vanishes 
otherwise. If the trispectrum is parity invariant (zero if l\ + h + h + h =odd), the cross power spectra C^ c = = C£ c 
vanish. The remaining power spectra and their covariance are described by the general forms of Eqn. (41), (42), and 
(44). 
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E. Gradient-Gradient 



The product of temperature gradients can in general be decomposed into three quadratic statistics 

[V,e a (n)] [V,e fc (n)] = t ab (h) 9lJ (h) + ]>>± i < b (n) ( 7±(n) +v(n)e ij (h) , (56) 

± 

where g^ is the metric on the 2-sphere, 

= ^ e T^)i(eeTeA , (57) 

gives the basis for a trace-free symmetric tensor field on the sky, and 

ey(n) = (e e )j{e 4> ) i - (e g ) l (e 4> ) j (58) 

gives the basis for a trace-free antisymmetric tensor field on the sky. The flat-sky versions of these statistics were first 
employed by [10] for CMB lcnsing and note that q, u, v are analogous to the similarly named Stokes parameters for 
polarization. 

As is the case for the CMB polarization, these three fields may be decomposed into multipole moments of the 
spherical harmonics and spin-2 spherical harmonics [9], 

^(n) = £if M F^(n), 

LM 

^(n) = £< M l^(n), 

LM 

[q ± iu) ab (n) = J> ± ib)f M ±2 F L M (n) , (59) 

LM 

where the moments follow the general prescription of Eqn. (38) with x = t, e, b, v and weights 

V hh (L) = l - y/h{h + \)y/h(h + l) ( !_\ \ o ) , Odd , 

ehh {L) = \ y/h(h + l) Vhih + l) ( \ ^ \ ) , even , 

b h i 2 {L) = - l -Vh(h + l)Vh(h + l) 1\ 2) ' ° dd ' (60) 

If the trispectrum is parity invariant, cross power spectra exist only among [t, e, g, s) and (f3,v,c). These power 
spectra and their covariance again are described by the general forms of Eqn. (41), (42), and (44). 



Ulio)' 6Ven 



F. Temperature-Hessian 

Similarly to the gradient-gradient case, the product of the temperature and the second derivatives or Hessian of 
the temperature field can be decomposed into three quadratic statistics 

e a (n)V 4 V,e b (n) = h ab (n) 9ij (n) + J> ± Mb]« 6 (n)<r±(n) , (61) 

± 

which themselves may be decomposed into multipole moments of the spherical harmonics and spin-2 spherical har- 
monics, 

h ab (h)=J2h a L b M Y L M (n), 

LM 

[ Vl + im] ab (n) = £(e + mfu ±2Y L M (n) , (62) 

LM 
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where the moments follow the general prescription of Eqn. (38) with x = h, e, f3 and weights 

h hh (L) = ~l 2 (l 2 + 1) ( q q q ^ , even, 
= ~k(l2 + l)SiiI 2 , 



eW2(L)= 2V(^2)! lo -2 2 J ' even ' 



^W^-W^^fn 1 'of] . odd, (63) 



2ty (/ 2 -2)! V° - 2 2 

Again parity invariance requires that power spectra exist only between (h, e, t, e, s) and ((3,b,v,c). Likewise the 
general formula for power spectra and their covariance again apply. 

G. Temperature- Temperature Hessian 

Auxiliary two-point statistics can be formed from the fundamental ones above. For example 

ViV^G^nje^n)] = t ab (h) 9ij (n) + ± i]°'(n)4(n) , (64) 



and 



t ab (n) = J2t a LMYl I (n), 

LM 

[q + l u] ab (n) = £(e + il)f M ±2 Y L M (h) , (65) 



LM 



imply that 



L(L + 1) (h l 2 L 




1 (L + 2)\ fh l 2 L 



^hh e hh e l r jh "1" e hh e Wi ^ y ^ — 2)! ^ J ' 

6iii 2 = ^1/2 + ^2/1 + Aii 2 + A 2 ii = • 

Again power spectra follow from the general relations. 

H. Temperature-Gradient Divergence 

The divergence of the temperature-gradient field of §IV D is also an auxiliary statistic 

V 4 [e a (n)V,e b (n)] = d ab (n)Yi 4 (n) , (66) 

with 

d Qb (n)=^< M F i M (n). (67) 

LM 

The weights are related to the others as 



d hh = y/L(L+lMh + l) (q 1 1\ 1 ) , even, 



= -iy/L(L + l)g hl3 = 2t hh + 2h hh . (68) 
Again power spectra follow from the general relations. 
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I. Cubic Statistics 

Finally the cross correlation of cubic statistics with linear statistics are also related to the quadratic statistics 
introduced above. For example 

(e^nOe^nOe 3 ^) 9 4 (n 2 )) = ]T(-im s(3) + ( Vr" l (Ai)*T(n 2 ) • (69) 

Im 

More generally, the cubic power spectra corresponding to the various x-statistics are given by 

C f {3) = E ^^^(^^(^(-^^^X 1 / 2 ^), (70) 

lll2l3L 

with Gaussian noise bias 

l± L 

and auxiliary noise variance 

^(3,12) = 1 yj^w [^sw + (-i^+'-^w] crcr , 

^3,34) = i. yj^tw + (-i) ii+ ' 2+L ^w] crcr* , (72) 

/ 1 L 

i.e. the multipolc index of the power spectra is no longer the diagonal of the trispectrum configuration but rather one 
of its sides. These cubic statistics thus probe a different projection of the trispectrum information but are based on 
the same set of filters and employ the same statistical formalism. 

V. CMB LENSING 

In this section, we consider the trispectrum signal generated by the weak gravitational lensing of the CMB tem- 
perature anisotropies by the large-scale structure in the Universe. In §VA we derive the full trispectrum for lensing 
and relate it to the underlying deflection (or convergence) power spectrum. Zaldarriaga [6] previously considered 
the lensing trispectrum in the small-scale or flat-sky approximation. The gravitational lensing effect is known to be 
dominated by potential fluctuations on the largest scales where an all-sky treatment of the trispectrum is desirable 
[7]. In §VB, we show that the signal-to-noisc in the trispectrum is both large and highly configuration dependent for 
experiments that can resolve multipole moments I > 1000. The divergence statistic introduced in §IVH is shown in 
§V C to be optimal for measuring the underlying deflection power spectrum at its peak at low multipoles. It benefits 
from substantially higher signal-to-noise as compared with the gradient-gradient quadratic statistics introduced by 
Zaldarriaga & Seljak [10]. Finally we consider tests for the robustness of the divergence statistic with differing filter 
sets in §VD and the degradation in power spectrum estimation from lensing covariance in §VE. 

A. General Trispectrum 

We begin by briefly reviewing the effect of gravitational lensing on the harmonics of the CMB temperature field and 
refer the reader to [7] for details of its calculation in a given cosmology. For reference, we employ the same ACDM 
cosmology used there, with parameters Q m — 0.35, Oa = 0.65, h — 0.65, n — 1 and Sh — 4.2 x 10~ 5 . 

Weak lensing of the CMB remaps the primary anisotropy according to the deflection angle V</>(n) 

6(n) = 9(n + V<f>) 

= e(n)+V^(n)V J 9(n) + ... , (73) 

where the tilde represents the unlcnscd field and . . . represent higher order terms in the Taylor expansion. The lensing 
potential field </>(n) is a lensing-probability weighted projection of the Newtonian potential along the line of sight [see 
[7], Eqn. (21)]. 
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FIG. 1. The power spectrum of the deflection angle in the fiducial ACDM model. Errors boxes represent the la errors from 
Gaussian noise on the divergence statistic binned in the bands shown. The divergence estimator of Eqn. (80)-(81) is optimal 
for the low multipoles and reduces the variance in the power spectrum estimation by more than an order of magnitude as 
compared with the gradient-gradient statistics of [10]. 



The spherical harmonic coefficients of the lenscd CMB temperature field become 
&im « e lm + J dnY7"*(n)V^(n)V*e(n) + ... 

= ®lm + EE <t>LMQl>m>{-lT ( ' , J^) F lLV +..., (74) 



LMl'm' 



where 



' " ;F lLl ' = I [L(L + 1) + I' (I' + 1) - 1(1 + 1)] ( ' ' L 



Y (2Z + 1)(2Z' + l)(2i + 1) 2 L v ; v ; v " \0 

= -VL(L + l)l'(l' + l) (J l \ , even, (75) 

where recall that "even" denotes the fact that only l + V + L=cvcn is non- vanishing. Gravitational lcnsing generates a 
change in the power spectrum that has been well studied [11-13,7]. It produces two changes to the 4 point function. 
The first is that the unlcnsed C\ in the Gaussian 4-point contribution must be replaced with the lensed C\. The 
second is that it generates a trispectrum an underlying reduced form [see Eqn. (23)] of 

T{£{L) = CfC l2 C u F h Li 2 F hLU . (76) 



Note the geometric interpretation: the lensing generates a trispectrum or quadrilateral configuration of l\ . . . U where 
one of the diagonals is supported by the lcnsing potential power spectrum Cff. Note that the power spectrum of 
the deflection field is given by L(L + l)Cff and is the fundamental quantity of interest. It is plotted in Fig. 1 for 
the fiducial ACDM cosmology. It is important to note that most of the power in the deflections is coming from 
a rather large scale or low multipole L w 50. Contrast this with the more familiar convergence power spectrum 

7Z» = [£(£ + i)/2] 2 cf 



C£ K = [L(L + l)/2] 2 Cf^ which peaks at much smaller angular scales. 



B. Total Signal-to-Noise 

From the considerations of §111 B, we can calculate the total signal-to- noise in the trispectrum for lensing in the 
full-sky formalism. Flat-sky estimates of the total signal-to- noise have been calculated in [6] (see also Appendix B). 
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500 1000 1500 2000 

k 

FIG. 2. Contributions to the (S/N) 2 from trispectra configurations with a fixed diagonal L and maximum side length ii, 
summed over the remaining three sides. Solid lines represent the full calculation of the trispectrum terms; dashed lines represent 
the pairwise approximation of Eqn. (77). The signal-to-noise in the low L trispectrum is highly dependent on the configuration. 

The all-sky expressions are cumbersome to calculate due to the presence of the Wigner-6j symbol that expresses 
the alternate recouplings of the trispectrum Vs. We use the recursion technique outlined in Appendix A for these 
calculations. 

In Figure 2, we show the signal-to-noise contributions in a given mode L = 50 of the trispectrum from a given l\ 
(summed over l 2 , h, h) assuming an ideal experiment Cf ot — C\. The signal-to- noise is quite high and approaches 
unity per l\ mode and L mode at h w 2000. Moreover, the contributions as a function of l\ show striking features. 
These features can be understood by approximating the trispectrum by its fundamental pairing (Zi, ^2)7 (^3^4) in 
Eqn. (23) 

T}£{L) « P}£ = Cf (c h F hLh + C h F hLh ) (c u F hLl4 + C h F ULh ) . (77) 

Figure 2 (dashed lines) verifies that this is a very good approximation for the range of interest h 3> L. The reason is 
that these configurations represent flattenned quadrilaterals where one diagonal is much greater than the other. Since 
lensing effects peak at low L, the other pairings are correspondingly suppressed. These properties are hidden in the 
real-space 4-point function and highlights the benefit of considering harmonic-space statistics. 

To the extent that Ci is constant, the two terms within each set of parenthesis in Eqn. (77) cancel. Thus, the 
trispectrum picks out features in the underlying unlensed power spectrum, specifically those due to the acoustic peaks 
in the power spectrum. Note that the effects on the power spectrum itself exhibits the same effect: lensing acts to 
smooth the acoustic features in the spectrum. This structure implies that optimizing the Z-range of filters in the 
quadratic statistics is important for maximizing the signal-to-noise. 

The cumulative signal-to-noise integrating out to l\ — Z max as a function of L in an ideal experiment is shown in 
Fig. 3. The approximation of Eqn. (77) begins to break down as L approaches Z max but always in the sense that it 
underestimates the total signal-to- noise (dashed lines vs. points). This breakdown occurs since the two diagonals of 
a trispectrum quadrilateral become comparable and either can be supported by the lensing power in Eqn. (76). We 
also show in Fig. 3 (solid lines) the cumulative signal-to-noisc for the Planck satellite [14] with Cf ot taken from [3]. 
Planck approximates an ideal experiment with Z max w 1600. 

Finally, under the approximation of Eqn. (77) which slightly underestimates the total signal-to-noise, we can plot 
the cumulative signal-to-noise summed over all L as a function of l max (see Fig. 4). Again an Z max ss 1600 approximates 
the Planck experiment whose total (S/N) 2 w 3100. 
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FIG. 3. Cumulative signal to noise in the trispectra configurations with the diagonal L summed over l\ . . . h. Dashed lines 
represent an ideal experiment where Ci = Cf ot out to a maximum I = i max ; solid lines represent the Planck experiment. Lines 
represent the approximation of Eqn. (77); points represent the calculation using the full trispectrum for an ideal experiment. 

C. Divergence Statistic 

The structure in the signal-to-noise curves imply that it is important to select a quadratic statistic that captures 
this structure. Following the considerations of §IVA, we can search for the quadratic statistic that optimizes the 
signal-to-noise ratio in the power spectrum of the deflection. Recall that in general the multipole moments of the 
quadratic statistic with filters /" and 



x?m = (-i) M E E <i 2 (L)eh mi ei 

Ixmi l 2 m 2 



2L + 1 h l 2 L 
" ""V ^ I mi m 2 —M 



(78) 



are defined in terms of the weight function a;"* . Under the approximation for the trispectrum of Eqn. (77), Eqn. (49) 
gives the optimal weights as 



12 C'i 2 Fi lL i 2 + C^F^Lh 

X hh ^ £ftot£<tot 
h h 



(79) 



Since the lensing trispectrum is symmetric in l\ <-> l 2 , the temperature-gradient divergence statistic whose weights are 



d hh - fl\fl 2 F h 



1 rl 



47T 



2L + 1 ' 

is optimal if the underlying temperature field is first filtered with 

A 



(80) 



ft = ft 
ft = ft 



(Jtot ' 

Ci 



We choose the proportionality constant to return the properly normalized deflection power spectrum 

-1 



A = ^L(L+l)(2L + 1) 



E 

hh 



(Ci 2 Fi lL i 2 + C^F^Lh) 

2^tot (J^ot 



(81) 



(82) 
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FIG. 4. Approximate total (S/N) 2 in the trispectrum for an ideal experiment out to / = i max and the Planck experiment. 
The Planck experiment approximates an ideal experiment of I « 1600 with a (S/N) 2 ~ 3100. 



The signal-to-noisc in this statistic is maximal under the assumption that the trispectrum can be approximated as 
Eqn. (77) as is the case for L < several hundred. The total signal-to-noise as calculated from Eqn. (45) for Planck is 
formally (S/N) 2 rts 4050. Compare this with the gradient-gradient statistic of [10]; for the unfiltered e ab (n) = £(n) 
the (S/N) 2 S3 135. 

The total signal-to-noise is allowed to exceed the maximum estimate of the previous section since the latter is 
strictly a lower limit. However the underlying approximation that the estimates at all L's are independent breaks 
down when integrating over a wide range in L causing a small reduction in the number of independent modes (see 
§IVA). 

The divergence statistic may also be used as a direct estimator of the deflection (or equivalently the convergence) 
field itself. Generalizing the argument of [10], one can think of the quadratic statistic x in Eqn. (78) as averaging 
over many independent (high-/ or small-scale) realizations of the unlenscd CMB field with a fixed realization of the 
large-scale deflection field 



( x °lm ) cmb = <PlmJ 47r (2L + l) 51 x hh {ph F hLh + Ci 2 F hU ^) , (L > 0) . (83) 

For the divergence statistic filtered as in Eqn. (81) all contributions add coherently so that, 

( d iV)cMB = QL+l) ^ 2C t l ot C t l ot ip hFl2hh + Cl 2 F hLl 2 ) 

= ^L(L+l)0 LM ■ (84) 

From the multipole moments one can reconstruct the deflection or convergence map [15]. Of course the fact that we 
average over only a finite number of independent modes of the primary anisotropy means that the resulting map will 
be noisy, with noise properties given by the Gaussian noise power spectra. 



D. Robustness Tests 



The divergence statistic contains enough signal-to-noisc for Planck that the data may be further subdivided to 
check for robustness of the statistic. Especially worrying is the possibility that galactic and extragalactic foregrounds 
and systematic effects might generate a false signal. Even if these contaminants contribute only Gaussian noise, one 
must subtract out the noise bias from the power spectrum estimators with the filters defined in Eqn. (81), 
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FIG. 5. Degradation in the total (S/N) 2 in the power spectrum due to covariance from gravitational lensing. The degradation 
is minimal for the Planck experiment or any that is cosmic variance limited only out to I ~ 2000. 

Recalling the discussion of the filters in §IVB, we can eliminate Gaussian noise bias as well as noise-correlation 
between differing L by defining non-overlapping filters sets (/*, f 2 ) and (/j 3 , ff). The resulting estimates of the 
deflection field d]^ M and d 3 ^ M would then have statistically independent Gaussian noise properties such that the noise 
bias is eliminated in their cross correlation. Furthermore, if the signal is really due to lensing the various estimates of 
the deflection power spectrum must agree within their errors. 

The price of dividing up the sample in this way is the signal-to-noise in any given set. For example, by band-limiting 
the filters of Eqn. (81) to 500 < Z < 1400 for the (12) set and I > 1400 for the (34) set, the total signal-to-noisc 
is reduced by ~ \pl and correspondingly the errors in Fig. 1 are increased by <~ Note that in this case, the 
underlying filtered temperature maps contain no power in the multipoles of interest L ~ 100. Such a scheme would 
still yield a highly significant detection and help protect against contaminants. Filtered versions of other quadratic 
statistics can also serve as consistency checks. For example, even the simple temperature-temperature s statistic of 
§IVC, yields a (S/N) 2 ~ 200 once it is is filtered according to the peaks in the Fig. 2. 

E. Power Spectrum Covariance 

One might worry that the high signal-to-noise in the trispectrum for lensing comes at the expense of degraded 
signal to noise in the power spectrum due to covariance between the estimators. Fortunately this is not the case in 
the I < 1000 regime of the acoustic peaks. From Eqn. (34)-(36), we can calculate the degradation in the total (S/N) 2 . 
This degradation is shown in Fig. 5 for the Planck satellite and an ideal experiment out to Z max - 

VI. DISCUSSION 

We have provided a systematic study of the angular trispectrum or four-point function of the CMB temperature 
field. Symmetry considerations dictate the fundamental form of the trispectrum and govern the Gaussian noise 
properties of its estimators. The large number of independent configurations of the trispectrum imply that even 
subtle physical effects may be detectable when all of the information in the trispectrum is brought to bear. 

In practice, extracting all of the information in the trispectrum will be a difficult computational task. We have 
thus also conducted a systematic study of the power spectra of quadratic statistics which probe different aspects of 
the trispectrum. Techniques developed for extracting power spectrum statistics from large data sets can then be 
brought to bear on the four-point function. The drawback is that this compression of information is in general a lossy 
procedure. We have therefore examined a wide range of quadratic statistics and prefiltering schemes. Given a target 
form for the trispectrum signal, these statistics can be optimized in their signal-to-noise for the power spectra. 



1G 



As an example, we have reconsidered the four-point correlation generated by weak lensing of the primary anisotropies 
by large-scale structure in the Universe as a means of recovering the power spectrum of the deflection angles or 
convergence [10]. CMB weak lensing provides a unique probe of the large-scale properties of the convergence field. 
We identify a specific quadratic statistic, the divergence of the temperature-weighted gradient field, that achieves the 
maximal signal-to-noise in this limit. For the Planck satellite, the total (S/N) 2 ss 4000 and represents a reduction in 
the noise variance on the convergence power spectrum by over an order of magnitude as compared with the gradient- 
gradient estimators of [10]. There is sufficient signal-to-noise to conduct filtering tests to eliminate noise-bias and 
check for consistency between multipole subsets of the data. We plan to explore further the properties of these 
estimators and their use in extracting cosmological parameters in a separate work. 

The lensing example illustrates the importance of examining the configuration properties of the trispectrum when 
designing statistical estimators based on the four-point function. Extracting the wealth of information potentially 
buried in the trispectrum will be a rich field for future studies. 

Acknowledgments: I thank A.R. Cooray and M. Zaldarriaga for many useful discussions. This work was supported 
by NASA ATP NAG5-10840 and a Sloan Foundation Fellowship. 

APPENDIX A: WIGNER-6J SYMBOL 
1. Useful Properties 

The Wigner-6j symbol expresses the relationship between two distinct couplings of three angular momenta 

Jd = Ja + Jb + Jc 
= Je J c 

= Ja + Jf, (Al) 

such that the eigenstates of the (ec) coupling are related to the eigenstates of the (a/) coupling as 

|(e C )d 7 > = |(a/)d 7 > v /(2 e + l)(2/ + l)(-l) s | I b d e f | , (A2) 

where Y. = a + b + c + d. Geometrically, the Wigner-6j represents a quadrilateral with sides (a,b,c,d) whose diagonals 
form the triangles (a,d,f), [b, c, /), (c,d,e), (a,b,e) or the three dimensional tetrahedron composed of these four 
triangles. It vanishes if the any of the triplets fail to satisfy the triangle rule. The symmetries are related to rotations 
of the tetrahedron that interchange the vertices. The result is that the symbol is invariant under the interchange of 
any two columns and under the interchange of the upper and lower arguments in any two columns. 
The Wigner-6j symbol can thus be used to permute the pairings in a set of Wigner-3j symbols 

Dv + D<-ir'— { : \) } _y (: j d _ s ) - (: *> i) (: ; d _ s ) . m 

or equivalently by the orthogonality relation of the Wigner-3j symbols 

{ : v, } - £i> i)e+ ' +,+ * (; 1 (, c i <A4) 

Finally, the Wigner-6j symbol obeys 

z> +1 >{"5/}{"5;}-5&t- (A5) 

and 
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2. Evaluation 



Closed form expressions exist for special cases of the arguments. For example 

a b e \ {-ly+b+e 



c d o J" ^TTP+T) ^' (A7) 

More generally, they may be computed efficiently by a recursive algorithm introduced by [16]. Let us define 

»&)-{***}. (A8) 

The h(ji) satisfy the recursion 

j 1 E( ]1 + l)h(h + 1) + F{ n )h( 3l ) + in + l)E{h)h{h - 1) = , (A9) 

where 

E(h) = {[.if (32 j 3 f] [(j2+.n + i) 2 -f 1 ] [jl-ik-hf] [{h + h + if-il}} 1 ' 2 
F{h) = (2ji + i){ji(ii + i)[-ii(ji + l) + h{h + i) +hih + i)] 

+h{h + l)b'i(ii + 1) +h(.h + 1) - .h(.h + 1)] 
+^ 3 (^3 + l)b'i(ii + l) - h(h + l)+h(h + l)] 

-2ji0'i + l)ii(ii + l)} (A10) 

For a stable recursion, one begins at both of the two ends ji m ; n = max(|j2 — J3I, \h — h\) Jimax = niin(j2 +.73, h + h) 
with the boundary conditions -E(jimin) = and E(ji max + 1) = and matches the two in the middle. 
The normalization is fixed by 



^(2 J1 + 1)(2? 1 + 1)K. ?1 )] 2 = 1, 

3l 

Sgn[Mjlma*)] = (-l)A+J3+^+'3 , (All) 

which follow from Eqn. (A5). 



APPENDIX B: FLAT SKY APPROXIMATION 



In the flat sky approximation, one decomposes the temperature field into Fourier harmonics 

(0(n 1 )...0(n n ))= J^-J ^<e(l 1 )...e(l n ))e il ^...e fl - fi '. (Bl) 

Statistical isotropy is enforced by demanding that the correlation function be invariant under an arbitrary translation 
and rotation in the plane. Parity invariance is enforced by demanding symmetry under inversion of the coordinates 
or reflection across one of the coordinate axes. 

As usual translational symmetry — > rL + C, where C is a constant vector, is enforced by the closure condition 
that the n-point function is proportional to 

(2tt) 2 (5(1 1 + . . . + 1„) . (B2) 

The wavevectors L thus form a geometric figure of n, possibly intersecting sides. Rotational invariance for the two- 
point function and rotational and parity invariance for the three-point function imply that the corresponding harmonic 
spectra are functions only of the lengths of the sides: 

(9(106(12)) = (27r) 2 ( 5(l 12 )C (il) , 
(6(106(12)6(13)) = (27r) 2 <5(l 12 3)i? (;i , i2 ,3) , (B3) 
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where = lj + . . . 4- lj and B should be symmetric against permutations of its arguments. For the four-point 
function, rotational and parity invariance implies that the quadrilateral formed by the four wavevectors is a function 
of the lengths of the sides and the lengths of the two diagonals, 

(9(1!) . . . 6(1 4 )) - {2irf [5(l 12 )d(h 4 )C h C l2 + 6(l 13 )6(h 4 )C h C l3 + 6(l 14 )8(h 3 )C h C u } 

H2n) 2 5(i 1234 )T^(i 12 ,i 13 ). (B4) 

To parallel our treatment of the all-sky four-point function let us break Q into its three distinct pairings and demand 
symmetry with respect to permutation of the arguments 

T S:;S = p tS^ + p (hi 3 ) + p lii:]^) . m 

Note that £14 is a function of the other two diagonals. P is symmetric under interchange of its upper and lower 
arguments as well as ordering within them. 

The relationship between the all sky and flat sky spectra can be obtained by noting that [7] 



y^]T^y™(n)e"^, (B6) 



where <j>\ is the polar angle of 1. It is a straightforward exercise to show that 



Ci = C, 



(0 ' 



Bi M3 = I o o o ) V ^ B (h,h,h) • (B7) 



For the trispectrum, we begin with the general correspondence 



(B8) 



where "c" denotes the subtraction of the Gaussian piece in Eqn. (B4). We then exploit the pair symmetry of the 
trispectrum exhibited in Eqn. (B5) by breaking the delta function into the corresponding pairs. For the (12), (34) 
pairing, 

<5(li234) = J d 2 L 6(1, +l 2 + L)S(l 3 + 1 4 — L) , (B9) 
such that L = l\2- Expanding the delta functions in spherical harmonics, 

1 21+1 



v ' m 1 ...m i LM \i=l,4 ' J 



47T 



x(-i) M (i; i? t) L* t* \a- ii 1 "" 



mi m,2 M J \m 3 m\ ~M 
Substituting back in and integrating over the polar angles, we obtain the general correspondence 

= 2 -^^k + l)...(2k + l)(^ l l L Q ) ('» h Q L Q ) pg$(L) , (Bll) 
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from which wc can construct the relation for T, \ 2 . 

L3L4 



We can now also make the correspondence between the signal-to-noise in the all-sky and flat-sky formalisms. The 
weighting of four-point terms that maximizes the signal to noise is [6] 

HV J - ] ' fd\ [ d% [ d% [ ePWx^) '^jaL' , (B12) 



NJ tot ir 24(2tt) 4 7 J J J V ' C\fC\fC\fC\f ' 

where / s k y is the fraction of sky covered by the sample. 

The square of T in this expression contains cross term in the P pairings. 

S\ 2 ( S\ 2 ( S\ 2 ( S\ 2 _ / S\ 2 / S \ 2 / S" x ~ 



Nj tot W (12,12) + U ) (13,13) + W (14,14) +2 UJ(1 2 ,1 3) +2 W (12,14) ^ U J (13,14) ' ' ^ ' 

The correspondence with the all-sky expression Eqn. (33) can be established by considering the terms pair-by-pair. 
For example for the (12,13) term, one expands the delta function in the (12), (34) pairing as above and inserts an 
additional delta function 

7rJ(li234) = 7T J d 2 L 13 S(h +l 3 - L 13 )5(h + h + Li 3 ) , (B14) 

with the understanding that 5(0) = V/(2tt) 2 = 1/tt. Expanding the delta functions in spherical harmonics we can 
integrate over azimuthal angles to obtain 

S \ _ /sky y y ^_-qMi2+Mi 3 ( h h L\2^\ ( h h L 12 



NJ no m 24 ^ ' ' \mi m 2 M 12 \m 3 m± -M 12 

( LZ < Ld ) Iimi...Um4 L 12 M 12 L 13 M 13 



h h L 13 \ ( l 2 h L 13 \ PiliT'iL^PililiLis) 



mi m 3 Mi 3 I \m 2 m 4 -M 13 J Cf ot Cf ot C} ot C} ot ' 

/ \ / fi Lo 13 U 



y h h h U 
jhh*{ t \ r>hh 



_ /sky \- V- (_ U h+l 3 i h h L 12 \ P lllT( L ^) P llll( L ^) /tdi r\ 

24 ^ ^ 1 J 1 k h L 13 J C, tot C, tot C, tot C, tot ' [ ' 

h...U L 12 L 13 ( h h h U 

where we have used Eqn. (A4) to rewrite the sum over the Wigner-3j symbols in terms of the 6j-symbol. Proceeding 
similarly for all terms in the signal-to-noise expression, we obtain 

S_\ 2 = /sky\p V- 1 \ T i 3 h( L )\ 2 / B1 „x 

Nj tnt 24 ^ ^ 2L+lC, tot C, tot C, tot C, tot ' 1 ; 

7 tot L hl 2 l 3 U ll ' 2 '3 '4 

where we have employed Eqn. (A6) to reexpress the (13, 14) term. The factor of 24 comes from the 4! permutations 
of each quadruplet in the all-sky expression. The / s k y term is the reduction in signal-to-noise due to incomplete sky 
coverage. 
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